Distinguishable RGE running effects between Dirac neutrinos and 
Majorana neutrinos with vanishing Majorana CP-violating phases 



Zhi-zhong Xing* and He Zhangt 

CCAST (World Laboratory), P.O. Box 8730, Beijing 100080, China 
and Institute of High Energy Physics, Chinese Academy of Sciences, 
P.O. Box 918, Beijing 100049, China 



Abstract 



In a novel parametrization of neutrino mixing and in the approximation of r- 
lepton dominance, we show that the one-loop renormalization-group equations 
(RGEs) of Dirac neutrinos are different from those of Majorana neutrinos 
even if two Majorana CP-violating phases vanish. As the latter can keep 
vanishing from the electroweak scale to the typical seesaw scale, it makes sense 
to distinguish between the RGE running effects of neutrino mixing parameters 
in Dirac and Majorana cases. The differences are found to be quite large in 
the minimal supersymmetric standard model with sizable tan (3, provided the 
masses of three neutrinos are nearly degenerate or have an inverted hierarchy. 
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I. INTRODUCTION 



Since 1998, a number of successful neutrino oscillation experiments [1-5] have provided 
us with very convincing evidence that neutrinos are massive and lepton flavors are mixed. 
A global fit of current experimental data [6] yields two neutrino mass-squared differences: 
Am 2 21 = m 2 -mj = (8.0 ± 0.3) x 10" 5 eV 2 and Am 2 2 = m 2 -m 2 = ±(2.5 ± 0.3) x 10~ 3 eV 2 . 
The upper bound of every neutrino mass is expected to be m { < 0.23 eV (for i = 1, 2, 3) [7]. 
Three neutrino mixing angles, associated respectively with solar, atmospheric and CHOOZ 
experiments, are well restricted: 30° < 9 12 < 38°, 36° < 9 23 < 54° and 9 13 < 10° at the 
99% confidence level [6]. But whether neutrinos are Dirac or Majorana particles remains 
an open question. If the neutrinoless double-beta (PP 0u ) decay is eventually observed, we 
shall make sure that neutrinos are Majorana particles. If there is no experimental signal 
for the j3j3 0u decay, however, we shall be unable to conclude that neutrinos are just Dirac 
particles [8]. Although most theorists believe that neutrinos should be Majorana fermions, 
there do exist some interesting models which treat massive neutrinos as Dirac fermions [9] . 
Before the nature of neutrinos is ultimately identified by the future neutrino experiments, 
it is worthwhile to study the phenomenology of both Dirac and Majorana neutrinos. 

Since many viable neutrino models are proposed at a superhigh energy scale, one has 
to take account of radiative corrections to their consequences on lepton flavor mixing and 
neutrino oscillations at low energy scales. Inversely, one may examine the evolution of those 
observed neutrino mixing parameters from the electroweak scale A EW ~ 10 2 GeV up to a 
superhigh energy scale by using their renormalization-group equations (RGEs). This kind 
of study is very useful, both for model building itself and for distinguishing different models. 
It is also useful to reveal the intrinsic differences between Dirac and Majorana neutrinos, as 
we shall explicitly show in this paper. 

A simple extension of the standard model (SM) or the minimal supersymmetric standard 
model (MSSM) is to introduce very heavy right-handed Majorana neutrinos and keep the 
Lagrangian of electroweak interactions invariant under SU(2) L xU(l) Y gauge transformation. 
The smallness of left-handed Majorana neutrino masses can then be explained via the well- 
known seesaw mechanism [10]. Below the typical seesaw scale A ss ~ 10 14 GeV, where heavy 
Majorana neutrinos become decoupled, the effective neutrino coupling matrix k obeys 



at the one-loop level [11], where t = ln(yu/Ass) with /i being an arbitrary renormalization 
scale between Aew an d Ass, Y x is the charged-lepton Yukawa coupling matrix, C = —1.5 



(SM) or C = 1 (MSSM), a M « -3g 2 + Qy 2 + A (SM) or a M « -1.2g 2 - Qg 2 + Qy 2 (MSSM). 



Here g 1 and g 2 are the gauge couplings, y t stands for the top-quark Yukawa coupling, and A 
denotes the Higgs self-coupling in the SM. 

If neutrinos are Dirac particles, their Yukawa coupling matrix Y v must be extremely 
suppressed in magnitude to reproduce the light neutrino masses of 0(1) eV or smaller at 
low energy scales. The running of Y u from a superhigh energy scale (e.g., A ss ) down to A EW 
is governed by the one-loop RGE [12] 




(1) 




(2) 
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where u = Y u Yj, a D « -0A5gf - 2.25g% + 3y 2 (SM) or a D « -0.6# 2 - 3#| + 3y 2 (MSSM). 
In writing out Eq. (2), we have safely neglected those tiny terms of 0(uj 2 ). 

One may use Eq. (1) or (2) to derive the explicit RGEs for neutrino mixing parameters 
in the flavor basis where Y l is diagonal and real (positive). In this basis, we have k = V M 7tV^ 
with k = Diagj/%, n 2 , ^3} for Majorana neutrinos; or cu = V D wV D with to = Diag{y 2 , y 2 , yf } 
for Dirac neutrinos. V M or Vd is just the lepton flavor mixing matrix. The unitarity violation 
of V M is extremely tiny and can safely be neglected in all realistic seesaw models [13]. At 
A EW , Majorana neutrino masses are m i = v 2 n i (SM) or m i = v 2 K i sin 2 j3 (MSSM) and Dirac 
neutrino masses are m i = vy i (SM) or m i = vy i sin (3 (MSSM) with v rs 174 GeV. A general 
parametrization of V M or Vd is V M = Q M UP M or V D = Q D UP D , where P M (or P D ) and Q M 
(or Q D ) are two diagonal phase matrices, and U is a unitary matrix containing three mixing 
angles and one CP-violating phase. Note that P D , Q B and Q M have no physical significance, 
but P M = Diag{e* p , e la , 1} consists of two non-trivial (physical) phases which are commonly 
referred to as the Majorana CP-violating phases. A novel parametrization of U is [14] 

/ Sis u c + c x c v e-^ s t c u c - c t s u e^ s t s \ 
U = c t s u c - Sjc^e - ** c { c v c + s x s v e-^ c x s , (3) 

V ~ S u S ~ C u S C J 

where c L = cos8 t , s v = sin^, c = cos^, and so on. This parametrization, together with the 
approximation of r-lepton dominance (i.e., YjY x Diag{0, 0, y 2 .} in view of y 2 <ti y 2 y 2 ), 
allows us to obtain the following RGEs for two Majorana phases p and a [12]: 



■ _ Cy 2 T r 
P 16n 2 



16tt 2 



7 2 2 ,7(22 2\ ,722 

Ci 2 c^s c (CT _ p) s ((T _ p) + C13 [s u s - c ) c p s p + ( 23 c u s c a s a \ , 

Cl 2 S^ 2 C ((J ^) S (--p) + ^3^ s2c p S P + C23( C ^ 2 - c2 ) C aS CT ] , ( 4 ) 



where p = dp/dt, a = dp/dt, = Ak^/ (k 2 — k 2 J, c a = cos a and s a = sin a (for a = p, a 
or a — p). A particularly interesting feature of Eq. (4) is that p = a = at a specific energy 
scale leads to p — a — 0, implying that p and a can keep vanishing at any energy scales 
between A EW and A ss . In this case, only three mixing angles (9 U 6 V) 6) and the so-called 
Dirac CP-violating phase undergo the RGE evolution. Note that a kind of underlying flavor 
symmetry may actually forbid two Majorana phases to take non-zero values in a concrete 
neutrino model. It is therefore meaningful to ask whether the RGE running behaviors of 
Majorana neutrinos with p = a = are identical to those of Dirac neutrinos. The purpose 
of this paper is just to answer such a question. 

In section II, we show that the one-loop RGEs of Majorana neutrinos with p = a = 
are analytically similar to those of Dirac neutrinos, but their expressions are not exactly 
identical. Section III is devoted to a numerical analysis of the RGE running behaviors of 
three mixing angles and the Dirac CP-violating phase, and to a careful comparison between 
the cases of Dirac and Majorana neutrinos. Four different neutrino mass spectra are taken 
into account in our calculations. We find that the differences between Dirac and Majorana 
neutrinos in their RGE running effects can be quite large in the MSSM with sizable tan /3, 
provided the masses of three neutrinos are nearly degenerate or have an inverted hierarchy. 
A brief summary of the main results is given in section IV. 
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II. RGES OF DIRAC AND MAJORANA NEUTRINOS 



The one-loop RGEs for three Yukawa coupling eigenvalues of Dirac neutrinos (y i with 
i = 1,2,3) and their four flavor mixing parameters {6 t , 9 V , 9 and 0) have been derived in 
Ref. [12]. Here we replace y { by m^. The RGEs of three neutrino masses, three mixing angles 
and one CP-violating phase can then be written as 



m 



ITLn 



16n 

m 



16tt 2 
16tt 2 



x - (a D + C y y u s 
(a D + Cy 2 T c 2 u s 2 ) , 
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where c, = cos and s. = sin 0. Note that the neutrino mass-squared differences m 



■mt = 



Am| x and m|- 



-m 



Am| 2 are much larger in magnitude than m 



•m; 



Am^, as indicated 



Am 2 x | 



|Am 2 2 | 



by current experimental data. Typically, Am 21 ~ 8.0 x 10~ 5 eV and 
2.5 x 10 -3 eV 2 [6]. Among three neutrino mixing angles, the RGE running of 9 V is expected 
to be most significant. The CP-violating phase may significantly evolve from one energy 
scale to another, if 9 t takes sufficiently small values. These qualitative features will become 
clearer in our subsequent numerical calculations. 

The one-loop RGEs for three effective coupling eigenvalues of Majorana neutrinos (/% 
with i — 1, 2, 3) and their six flavor mixing parameters (9 t , 9 U , 9, 0, p and a) can be found in 
Ref. [12]. Here we replace K t by m i and take p = a = at either A EW or A ss . As pointed out 
in section I, Eq. (4) assures two Majorana phase p and a to keep vanishing at any energy 
scales between A EW and A ss . One may safely simplify the RGEs of 9 t , 9 U , 9 and obtained 
in Ref. [12] for Majorana neutrinos by setting p = o = 0, and then compare them with their 
Dirac counterparts on the same footing. In this case, we arrive at 
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(m 3 - m 1 ) (m 3 - ,,c 2 

As a consequence of Am 21 <C |Am 31 | ~ |Am 32 |, the mixing angle Q v is most sensitive to 
radiative corrections. The RGE evolution of the CP-violating phase depends strongly on 
the smallness of Q h on the other hand. These qualitative features are essentially analogous 
to what we have pointed out for Dirac neutrinos. 

It is interesting to note that Eq. (7) can actually be obtained from Eq. (5) with the 
replacements a D a M and C =^ 2C, while Eq. (8) can be achieved from Eq. (6) with 
the replacements mf =>- m i (for i = 1,2,3). These similarities and differences imply that 
it is very non-trivial to distinguish between the RGE running behaviors of Dirac neutrinos 
and Majorana neutrinos with vanishing Majorana CP-violating phases. 

A rephasing-invariant description of leptonic CP violation in neutrino oscillations is to 
make use of the Jarlskog parameter J [15]. It explicitly reads as J = c l s l c u s v cs 2 s <t> in the 
parametrization taken in Eq. (3) for either Dirac or Majorana neutrinos. The one-loop RGE 
of J has been derived in Ref. [12]. Taking p = a = in the Majorana case, we obtain a 
simplified expression of j M , 

Cyl 
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167T 2 
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which is very analogous to j B of Dirac neutrinos, 
Cyl 



+ c 
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16vr 2 



( 2 _ 2 
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S n i 



l 2 — TTl^ - ' /#i< 3 — nt-^ 

It is obvious that Eq. (10) can be obtained from Eq. (9) with the replacements =^ m 2 
(for i = 1,2,3). Note that j D oc J D (or j M oc J M ) holds. This result implies that the 
Jarlskog parameter will keep vanishing at any energy scales between A EW and A ss , if it 
initially vanishes at either A EW or A ss . 

Two comments are in order. First, three mixing angles 9 t , 6 V and 6 in our parametrization 
can simply be related to 9 12 , 9 23 and 9 13 in the standard parametrization advocated by the 
Particle Data Group [16]. The relations are [17] 
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(9) 



(10) 



012 ~ 0; 



6 23 6 







13 



9 lS m9 (11) 

in the leading-order approximation. Thus 30° <0 V < 38°, 36° < Q < 54° and 9 l < 17° are 
expected to hold, as indicated by the global fit done in Ref. [6]. Second, the approximation 
of r-lepton dominance used in our analytical calculations is very reliable. In other words, 
the contributions of y 2 e and y 2 to all of the RGEs listed above are negligibly small. This 
observation has been verified by our numerical calculations, in which there is no special 
assumption or approximation. 
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III. NUMERICAL ILLUSTRATION AND COMPARISON 



In view of the fact that the absolute mass scale of three light neutrinos and the sign of 
Am\ 2 remain unknown at present, let us consider four typical patterns of the neutrino mass 
spectrum: 

• Normal hierarchy (NH): m 1 <C m 2 <C m 3 . For simplicity, we typically take m 1 = at 
A EW in our numerical calculations. Then m 2 = ^Am^ and m 3 = Am\ 2 \ + Am 21 
can be determined from current experimental data. 

• Inverted hierarchy (IH): m 3 <C m 1 < m 2 . For simplicity, w typically take m 3 = at 
A EW in our numerical calculations. Then m 2 = ^J\Am 32 \ and m 1 = ^J\Am 32 \ — Am 21 
can be determined from current experimental data. 

• Near degeneracy (ND) with Ara§ 2 > 0: m 1 < m 2 < m 3 . For simplicity, we typically 
take m 1 = 0.2 eV at A EW in our numerical calculations. 

• Near degeneracy (ND) with Am 32 < 0: m 3 < m 1 ^ m 2 . For simplicity, we typically 
take m 1 = 0.2 eV at A EW in our numerical calculations. 

In addition, we take Am\ x = 8.0 x 10~ 5 eV 2 , \Am 2 32 \ = 2.5 x 10~ 3 eV 2 , V = 33.8°, 9 = 45°, 
6 l = 0.5° and </> = 90° as typical inputs at A EW ~ M z in our numerical calculations. 

A. Neutrino masses 

In either the SM or the MSSM with small tan/5, the RGE running behaviors of three 
neutrino masses are dominated by a D or a M . The y 2 -associated term of m 2 in Eq. (5) or 
(7) becomes important only when tan/3 takes sufficiently large values in the MSSM [18]. 
Note that a M = 2a D holds in the MSSM, in which the running effects of m i for Majorana 
neutrinos are twice as large as those for Dirac neutrinos. 

The first plot in Fig. 1 illustrates the ratios Ri = m^ji) / m^M z ) changing with the 
energy scale fi in the SM for Dirac and Majorana neutrinos, where m 1 (M z ) = 0.2 eV and 
M H = 180 GeV (the Higgs mass) have typically been input. Since the running of m i is 
governed by a D or a M , R x ~ R 2 ~ R 3 holds to a high degree of accuracy. Furthermore, the 
behaviors of R i are actually independent of the initial value of m l and possible patterns of 
the neutrino mass spectrum. We observe that i?j in the Majorana case is always larger than 
R i in the Dirac case, and their discrepancy can be as large as 0.7 at fi ~ A ss ~ 10 14 GeV. 

The relation R x m R 2 « R 3 is also a very good approximation in the MSSM with small 
tan/3, as shown by the second plot in Fig. 1, where tan/3 = 10 has been input. It is clear 
that R { in the Dirac case is numerically distinguishable from R i in the Majorana case, in 
particular when the energy scale \x far exceeds M z . 

If tan P is sufficiently large, the common scaling of three neutrino masses in the RGE 
evolution will fail [18]. The splitting of R l , R 2 and R 3 , which increases with the energy scale 
fi, is illustrated by the third plot in Fig. 1 with the input tan/3 = 50. One can see that R i 
in the Dirac case is always smaller than R i in the Majorana case, and their discrepancy is 
distinguishable at the scales ji 3> M z . 
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B. Neutrino mixing parameters 



Radiative corrections to three neutrino mixing angles, the Dirac CP-violating phase and 
the Jarlskog parameter are all controlled by the r-lepton Yukawa coupling eigenvalue y T . 
Because of y 2 /(8n 2 ) « 1.3 x 1(T 6 (SM) or y 2 /(8n 2 ) « 1.3 x 1(T 6 (1 + tan 2 [3) (MSSM) 
at M z , significant RGE running effects are expected to appear in the MSSM case when 
tan/3 is sufficiently large [19]. To illustrate, here we simply concentrate on the MSSM 
with tan/3 = 50 and consider four typical patterns of the neutrino mass spectrum in our 
subsequent discussions and calculations. 

(1) In the NH case with m 1 = 0, the RGEs of 9 h 9 U , 9, and J can be simplified as 

a , C Vr a Cy 2 T ( 2 2 \ 

1 = + 16V C - S - CC 't> r ' 6 " = ~I6^ C " S " \ S ~ ° r ) ' 

Cy r (1 , \ 1 (2 2\ -1 -1 

J = ilv J I 2 ( sy - c ') " ( c2 " c ' s2 ) r l ' (12) 



in which r = 2m 2 /(m| — m 2 ) for Dirac neutrinos or r = 2m 2 /(m 3 — m 2 ) for Majorana 
neutrinos. Current experimental data yield r D rs 0.06 and r M ps 0.4. Both of them are too 
small to compensate for the strong suppression induced by y 2 T in Eq. (12). Thus the RGE 
corrections to those flavor mixing and CP-violating parameters are not important in the 
NH case. Note, however, that the radiative correction to can be very significant when 9 l 
is extremely small or becomes vanishing. We find that </> quickly approaches its quasi-fixed 
point 0q F = or 7r in the 9 l — > limit, an interesting phenomenon which is remarkably 
different from the non-trivial quasi-fixed point of discovered in the general (p 7^ o 7^ 0) 
case for Majorana neutrinos [20]. One can also see that both J = and J = hold when 
9 X vanishes; i.e., CP is a good symmetry in this limit. 

(2) In the IH case with m 3 = 0, we arrive at 

^ = = 0, 9 V = c v s v s 2 r' , 9 = +^ cs , 

j = [3c 2 - 1 + (si - cl) s 2 r'] , (13) 

where r' = (m| + m 2 )/(m 2 — m\) for Dirac neutrinos or r' = (m 2 + m 1 )/(m 2 — mj for 
Majorana neutrinos. We observe that radiative corrections to 9 { and are vanishingly 
small, and the correction to 9 is also insignificant. Nevertheless, the RGE running effects 
of 9 V and J may get enhanced by r', whose typical value reads r' D pa 60 or r' M ~ 120 at 
M z . Fig. 2 illustrates the evolution of 9 V and J in the IH case. The discrepancy between 
Dirac and Majorana cases is obviously distinguishable for both parameters, when the energy 
scale is much higher than M z . In particular, J7" D ~ 2J M holds at /i ~ A ss , because the 
corresponding value of 9 V for Majorana neutrinos is roughly half of that for Dirac neutrinos. 

(3) In the ND case with Am| 2 > and m 1 = 0.2 eV, the RGE corrections to those 
neutrino mixing parameters can significantly be enhanced by the ratios (mf + m 2 ) / (m 2 —m 2 ) 
in Eqs. (6) and (10) for Dirac neutrinos, or by the ratios (m^ + m-)/{m i — m-) in Eqs. (8) 
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and (9) for Majorana neutrinos. We illustrate the typical evolution behaviors of 9 h 9 U , 9 
and in Fig. 3. One can see that Majorana neutrinos undergo the RGE corrections more 
significantly than Dirac neutrinos. The discrepancy between these two cases is about 10° 
for either 9 or at /i 3> M z . It is therefore possible to distinguish the running of Majorana 
neutrinos from that of Dirac neutrinos. The difference between J D and J M is insignificant 
even at p ~ A ss , as shown in Fig. 4, partly because the increase (or decrease) of 9 l can 
somehow compensate for the decrease (or increase) of 9 and in the Majorana (or Dirac) 
case. 

(4) In the ND case with Ar/i^ < and = 0.2 eV, we get similar enhancements in 
the RGEs of those neutrino mixing parameters induced by the ratios (mf + m 2 j )/(mf — mf) 
for Dirac neutrinos, or by (m i + mf)j{m i — mf) for Majorana neutrinos. However, only the 
running of 9 is sensitive to the sign flip of Am| 2 , as one can see from Eqs. (6) and (8)-(10), 
in which 9 V and J are dominated by the term proportional to (m\ + mf) / {m\ — mf ) (Dirac) 
or (m 2 + m 1 )/{m 2 — m 1 ) (Majorana). Then the numerical results for 9 t , 9 U , and J in the 
present case are very similar to those in the ND case with Am^ > 0. For simplicity, we only 
illustrate the evolution of 9 in Fig. 5 by taking Am| 2 < 0. It is obvious that the running 
behavior of 9 for either Dirac or Majorana neutrinos in Fig. 5 is essentially opposite (or 
complementary) to that in Fig. 3, just due to the sign flip of Am§ 2 . 

IV. SUMMARY 

The main goal of this paper is to examine whether the RGE running behaviors of Majo- 
rana neutrinos are still different from those of Dirac neutrinos, if two Majorana CP- violating 
phases vanish at a given energy scale. For this purpose, it is essential to choose a suitable 
parametrization of the 3x3 lepton flavor mixing matrix, such that its two Majorana phases 
keep vanishing in the RGE evolution from one scale to another. We have pointed out that 
the novel parametrization used in Ref. [12], which consists of the mixing angles {9 l1 9 u ,9) 
and the CP- violating phases (0, p, a), does fulfill this requirement. Taking p = a = at the 
electroweak scale, we have carefully compared the similarities and differences between the 
RGEs of 9[, 9 U , 9 and for Majorana neutrinos and those for Dirac neutrinos. Our numerical 
calculations show that it is possible to distinguish between these two cases in the MSSM 
with sizable tan (3, in particular when the masses of three neutrinos are nearly degenerate 
or have an inverted hierarchy. 

Of course, the numerical examples presented in this work are mainly for the purpose 
of illustration. The point is that the nature of neutrinos determines their RGE running 
behaviors, and the latter may be crucial for building a realistic neutrino model. We expect 
that our analysis can not only complement those previous studies of radiative corrections to 
the physical parameters of Dirac and Majorana neutrinos, but also help us understand the 
dynamical role of Majorana phases in a more general picture of flavor physics. 
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FIG. 1. The running neutrino mass ratios R { = m^^/m^M^ (for i = 1,2,3), where the 
dashed and solid curves stand respectively for the Dirac and Majorana cases. 
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FIG. 2. The running behaviors of 9 V and J in the IH case with tan/3 = 50 and m 3 = at 
M z within the MSSM, where the dashed and solid curves stand respectively for the Dirac and 
Majorana cases, and J D {M Z ) = J M {M Z ) » 0.0014. 




FIG. 3. The running behaviors of 9 t , U , 9 and <f> in the ND case with Am| > 0, tan/3 = 50 
and m x {M z ) = 0.2 eV within the MSSM, where the dashed and solid curves stand respectively for 
the Dirac and Majorana cases. 
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FIG. 4. The running behavior of J in the ND case with Am^ > 0, tan/? = 50 and 
m 1 (M z ) = 0.2 eV within the MSSM, where the dashed and solid curves stand respectively for 
the Dirac and Majorana cases, and J D (M Z ) = J M (M Z ) « 0.0014. 




FIG. 5. The running behavior of 9 in the ND case with Am^ < 0, tan/3 = 50 and 
m 1 (M z ) = 0.2 eV within the MSSM, where the dashed and solid curves stand respectively for 
the Dirac and Majorana cases. 
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